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Abstract. There have been derived equations describing the static distributions of temperature and wind velocity at the 
transition region within the framework of the magnetohydrodynamics (MHD) of fully ionized hydrogen plasma . We have also 
calculated the width of the transition between the chromosphere and corona as a self-induced opacity of the high-frequency 
Alfven waves (AWs). The domain wall is a direct consequence of the self-consistent MHD treatment of AWs propagation. We 
predict considerable spectral density of the high-frequency AWs in the photosphere. The idea that Alfven waves might heat 
fV^ the solar corona belong to Alfven — we simply derived the corresponding MHD equations. The comparison of the solutions 

rf^ • to those equations with the observational/measured data will be crucial for revealing the heating mechanism. The analysis of 

those solutions will explain how Alfven waves brick unto the corona and dissipate their energy there. 
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1. ALFVEN MODEL FOR CORONAL HEATING 

The discovery of the lines of the multiply ionized iron in the solar corona spectrum [ lj posed an important problem for 
the fundamental physics — what is the mechanism of heating the solar corona and why the temperature of the corona 
is 100 times larger then the temperature of the photosphere. 

The first idea by Alfven J2l was that Alfven waves (AWs) JH] are the mechanism for heating the corona. AWs 
are generated by the turbulence in the convection zone and propagate along the magnetic field lines. Absorption is 
| proportional to a> 2 and the heating comes from the high-frequency AWs. The low frequency AWs reach the Earth orbit 

and thanks to the magnetometers on the various satellites we "hear" the basses of the great symphony of the solar 
turbulence. 

The purpose of the present work is to examine whether the initial idea is correct and to derive the MHD equations 
which give the dependence of the temperature on the height T(x) and the related velocity of the solar wind U(x). 
Our starting point are the MHD equations for the velocity field v and magnetic field B 



<? r p+divj = 0, j = pv, (1) 
sJ ' ( 1 B 2 \ 

c3: c, '( v 2 py2+e+ 2^) +divq=0 ' (2) 

5,(pv) + V-n = 0, (3) 



where 



q = p Qv 2 +/ijv + v-n visc ->fVr + S (4) 

is the energy density flux, p is the mass density, e is the internal energy density, x. is the thermal conductivity, h is the 
enthalpy per unit mass; 

S = — [B x (v x B) - v m B x (V x B)] , (5) 
Mo 
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is the Pointing vector and v m = c 2 £oPn is the magnetic diffusion determined by the Ohmic resistance pn and vacuum 
susceptibility £o; vacuum permeability is For hot enough plasma pn is negligible and we ignore it hereafter. The 
total momentum flux 

n = p w + p + n visc + n Maxw (6) 

is a sum of the inviscid part p vv + P of the fluid, with pressure tensor P = p\, 

nr = - 77 (dm + d kVi - 2 - s ik v ■ v) - Cfev -v, (?) 

the viscous part of the stress tensor, with viscosity 77 and second viscosity and lastly, the Maxwell stress tensor 

_ n Maxw = _^/ BA _l B 2 5 ^ j (8) 

with 8% the Kronecker delta. We model coronal plasma with completely ionized hydrogen plasma 

T 5/2 mJ 2 T 5 l 2 „ 

x = 0. 9t — 77=0.4^-^, C«0, (9) 

fr D T\ 1 4ne 2 



A = ln 



where ^ e is the electron charge, m e is the mass of electron, m ? is the proton mass, T is the temperature and « to t = n e + n p 
is the total number density of electrons and protons; p = m p « p . We suppose that /io = 4-n and £o = 1 /47T, but in the 
practical system all formulas are the same; as well as in Heaviside-Lorentz units where Hq = 1 and £0 = 1. As we 
mentioned above 

c 2 e 2 m l J 2 A 77 OAT 5 / 2 

^s^^rs?^ (10) 

i.e., the hot hydrogen plasma is sticky, dilute, and "superconducting," v m w 0. Let us mention also the relations 
xpa = l.5T/ql and 77/^ ps \^m^m v , 

1 e 2 my 2 A 

2. MAGNETOHYDRODYNAMIC EQUATIONS AND ENERGY FLUXES 

The time derivative <3,B which implicitly participate in the energy conservation Eq. (f2|i at zero Ohmic resistivity obeys 
the equation 

d,B = B-Vv, d, = <9,+v-V. (12) 
Analogously the momentum equation ([3]l can be rewritten by the substantial derivative 

p d,v; = -dtp + d k |tj (d k Vi + djv k - ^8 ik djv^j | + d t (£djVj) -i(Bx rotB),- . (13) 

In our model we consider AWs propagating along magnetic field lines Bo. We focus our attention on the narrow 
transition layer, where the static magnetic field is almost homogeneous and the waves are within acceptable accuracy 
one-dimensional. For the velocity and magnetic fields we assume 

y(t,x) = U(x)x + u(t,x)z, 

B(t,x) = B \ + b(t,x)z, (14) 

with homogeneous magnetic field Bqx perpendicular to the surface of the Sun. The transverse wave amplitudes of the 
velocity u(t,x) and magnetic field b(t,x) we represent with the Fourier integrals 

u(t,x) = — f u{co,x)e- im 'dco, (15) 



27T. 

2kJ- 



t,x) = — I b((0,x)e- imt d(0. (16) 



For illustrative purposes it is convenient to consider monochromatic AWs with u(t,x) = «(x)e imt and b(t,x) = 
b{x)t- imt . 

2.1. Wave equations 

For linearized waves the general MHD equations (fT3l l and (fl2t give the following system for u(x) and b{x) = 
b(x)/B 

(-i(0 + Ud x ) a = Vld x b+ -d x (r]d x a) , (17) 

P 

-icot = d x u-d x (ut), (18) 

where 

Va(x)=Bo/VvoP(x) (19) 
is the Alfven velocity. In our numerical analysis we solve the first order linear set of equations 

-id*|¥> = -^-lv1|¥> = K|¥}, (20) 



where w = d x u, and 




(VI = 1^= [u*,b ,w" 



/ -v k U 

M= v k (-ico+d x f/) -v k 

\ iaU -Vfc-im+dJJ) (vl-U 2 ) + ^d x 7] 

For homogeneous medium with constant 17, p, V"a, and £/, in short for constant wavevector matrix K, the exponential 
substitution V <x exp(ifct) in Eq. ( f20l > or equivalently Eqs. (fTTT i and (fT8l gives the secular equation 

iv k I/det(K-fel) = (Od (coq + iv k /t 2 ) - V^k 2 = 0, (21) 



where ©d = o — ^ is the Doppler shifted wave frequency. This secular equation gives the well-known dispersion 

Vlk 



Od (fi>D + iVk^ 2 ) = °f me AWs. This equation is quadratic with respect to 0) and cubic with respect to k. 



2.2. Wind variables 

We solve the wave equation ( f20b from "Sun' surface" x = to some distance large enough x — I, where the short- 
wavelength AWs are almost absorbed. This distance is much bigger than the width of the transition layer X, but much 
smaller than the solar radius. The considered one-dimensional < x < I time-independent problem has three integrals 
corresponding to the three conservation laws related to mass, energy and momentum. The mass conservation law 
Eq. (HJ gives the constant flow 

j = p (x)U(x) = PqUq = piUi = const, (22) 

where po = p (0) , p/ = p (I) , Uo — U (0) , and U\ =U(l). The energy conservation law reduces to a constant flux along 
the x-axis 

9 = <7wmM + fa) =pu(^U 2 + h\+q= const. (23) 

Here the first term describes the energy of the ideal wind, i.e., a wind from an ideal (inviscid) fluid. The second term 
q(x) includes all other energy fluxes; in our notations tilde will denote sum of the non-ideal (dissipative) terms of the 



wind and wave terms. In detail the non-ideal part of the energy flux q(x) consists of: the wave kinetic energy «|#| 
viscosity (wind °c| rj + £ and wave 0=77 components), heat conductivity and Pointing vector °S* t 



q(x) 



J 1 ~i2 

4M ■ 



%Ud x U 



\r}d x \u\ 2 - 



xd x T 



1 

2^0 



U\b\ -B Q Re(b*u) 



(24) 



where £ = i 7] + £. Here time averaged energy flux is represented by the amplitudes of the monochromatic oscillations, 

'(Reuf 



this is a standard procedure for alternating current processes. In our case we have, for example, (u 2 ) t - 

( ^ (u + u*) 2 ) t = \ \u\ 2 . The other terms from Eq. (0) are averaged in a similar way in the equation above 
The momentum conservation law Eq. © gives constant flux II = n vv 



n = i^ d (x)+u(x) = P uu+p- 

the sum of the ideal wind fluid and the other terms 

1 ,<m2 



n, 



II(X)EE 



4^0 



b -& X U, 



(25) 



(26) 



which take into account the wave part of the Maxwell stress tensor <^b 2 and viscosity of the wind . 

We have to solve the hydrodynamic problem for calculation of wind velocity and temperature at known energy and 
momentum fluxes. The problem is formally reduced to analogous one for a jet engine, cf. Ref. |4J] . We approximate the 
corona as fully ionized hydrogen plasma, i.e., electrically neutral mixture of electrons and protons. The experimental 
data tell us that proton temperature T p is higher than electron one T e . This is an important hint that heating goes trough 
the viscosity determined mainly by protons. However, for illustration purpose and simplicity we assume proton and 

T V = T. For such an ideal (in thermodynamic sense) gas the local sound velocity 



electron temperatures to be equal T e 

is 



Cp_ p 

c v p 



c 2 s (x) 
(m) = ^ 
p = n tot T = 



PT 



n) 

1 

U (m) 



IE. 

cv 



(27) 



p- 



'p " 
T 



1 

2»tot, 



where, as we mentioned earlier, h is the enthalpy per unit mass and e is the density of the internal energy. 

In order to alleviate the final formulas we introduce two dimensionless variables a and T which represent the non- 
ideal part of the energy and momentum flux respectively 



_ q(x) 



Po^o 



t{x) = 



Po^o 



and analogously for the wind velocity and temperature 

©(*) 



m- uix) 



T(x) 



U 



(m) U, 



2 ' 



©o = ®(0), 



(28) 



(29) 



where Uq = U (0). The energy and momentum constant fluxes Eq. d23T l and Eq. (l25T l in the new notation take the form 



Pot/ 3 

n - n(o) 

PoU 2 



- -zU +c p @- 



c p @o, 



U + &/U-T=l+& . 



(30) 
(31) 



From the second equation we express the dimensionless temperature and substitute in the first one. Solving the 
quadratic equation for the wind velocity U we derive 



U(x) = — j-y (j+s 2 + yT(x) - y/@{x)j 



(32) 



where for the discriminant we have 



{s 2 - l) 2 - 2(y2 - l)a 2 + jt [jt + 2(1 +s 2 )] , 
2_C 2 (0) _ 2 , . (dp\ _ yr(x) 



Here 7 is the constant ratio of the heat capacities, and s = c s (0) /Uq is the ratio of the sound and wind velocity at x = 0. 
We suppose that initial wind velocity is very small 1/(0) <§C c s (0). The velocity distribution Eq. (l32l can be substituted 
in Eq. OTT ) and we derive a dimensionless equation for the temperature distribution 

&(x) = U{x) (\ + s — + yr(jc) - E/(jc)^ , (34) 

The solutions for velocity U(x) Eq. d32b and temperature T(x) Eq. ( [34-b distributions are important ingredients in our 
analysis and derivation of a selfconsistent picture of the solar wind. 



2.3. Entropy generation and equation for the temperature distribution 

Dissipated power density is equal to p T times the substantial derivative of entropy per unit mass 

pTd t Sf = div(>cVT)+pn(—) +o' lk d kVi . (35) 

From the right-hand side of this equation for the heat transfer the first term describes the thermal conductivity, the 
second one is the Ohmic power which is negligible pa ~ 0, and the last term gives the heating by viscous part of the 
strength tensor 

a ik = n ( d i v k + d k Vi - l^divv ) + C<5;A-divv. (36) 
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The time averaged density of viscous heating power is a sum of wind and wave parts 

{al k d k Vi) t = « (d x Uf + l - n \M 2 ^^ (lW 2 + \rf |#| 2 ) , (37) 

where 

— 4_ - 

<i; = -77 + C with pn = 0. 

We solve the transition layer problem in a space interval < x < I and introduce convenient for further work 
dimensionless variables 

_ x — — d x U d x u 

x = - , W = djU = , w = drU = , 

/ IU Q ' IU 

it 7 C _ (m) U$x 

77 = , l = , x= 5—. (3a) 

' po/f/o S Po/f/o lp U$ 

Analogously introducing reduced temperature and dimensionless temperature gradient F 

- T © _ _ 

T = — = — , F = d x T, F = djT. (39) 
we rewrite the heating power density by heat conductivity as 

div(^Vr) = PqU ° &0 [(djK)F + xdjF] . (40) 



The hot corona we model as ideal gas 



V 

S = N In — + Nc v In T + N {cy + const) , 



(41) 



of electrons and protons 



The substitution here n e =n 



Se ~t~ Sp 



m p Np + m e N e 



N p = n p V, N e = n e V. 



p — | n tot gives for the time derivative 

(d ( y(4=T-vy = ^ 

(m) 



d^tot 

n tot 



-c v - 



(42) 



where time averaging is over wave period. In such a way taking into account the mass conservation law j = (m) n tot U 
we obtain for the irreversible heating power 



P Td t y = —( P u = j) 

(m) 



d x U 



3\ d x T 



CV -2J T 



p t/ 3 ®o fW , c v F 



I 



U 



(43) 



The substitution of this equation together with Eq. d40b and Eq. ( 1371 ) in Eq. ( |3"5T l gives an explicit expression dimen- 
sionless temperature gradient 



_— _ — WT C — 2 tT ,^|2 

^F = (c v -d,x)F + - r --W - Wq \w\. 

The thermal conductivity Eq. © using the dimensionless variables can be rewritten as 

Ao ^5/2 _ J7o 



no 



PolUo' 



A simple differentiation gives 



_ _ Ao-5/2 _ (m)xb 
x=x —T', xr = — — , 
A PolUo 

3 _ 1 _ , 
A = Ao + -lnr+-lnC7, A =A(x = 0). 



3F 1W 
d Y A = -= + - = , 
2 T 2U 



(44) 

(45) 
(46) 



dx* _ d^7 
H rf 



5 3 \ F 1 W 



(47) 
(48) 



Then the gradient of the kinetic coefficient in Eq. (l44l can be expressed by its logarithmic derivative and we have not 
differentiations in the equation which have to be solved numerically. For the wind velocity we can proceed in a similar 
manner. 



2.4. Momentum equation for the wind velocity 

The x-component of momentum equations (0 and ( [TBI can be also rewritten as 



p(d,U + Ud x U) = -d x 



d x ($d x U), 



(49) 



where we have gradient from the sum of pressure and time averaged magnetic pressure of the propagating AWs. Now 
we have to express all variables by dimensionless ones 



p = n tot T = p Q Uo&o=, U = U U, b = B b, 



where 



x = lx, %=p lU %, V A = U V A , P 

_ Va(0) 



V a = 



o 



= alio V U, 



/MOP U 
is the Alfven velocity. Using the relation which can be easily checked 

tt2 



= £ = = L 

Po U' 

= V A (0) 



~ p U 2 U(x) 1/(0) 



the momentum equation takes the form 

|dyW = (l -d^ ) 
Then from the second row of Eq. (l20l we have 



-\- , f®oT a 2 t 2 
^ + — b 
U 4 



Ud x b= (i03-W)b + w, W = d x U, 



and easily derive 



Now we can substitute dv 



Ud-, 



= -2W 



- b w + bw* . 



in the momentum equation d52l i which finally takes the derivatives free form 



(TW 




a I — 




2 ^* , 




-S) 


+ -= [2W 


b 


— bw 






4U V 







(50) 



(51) 



(52) 



(53) 



(54) 



(55) 



where we taken into account that U(0) = 1. Let us mention also the relation between plasma j3 -parameter and 
dimensionless variables 



Bl/2no r "[/' ^ 7a 2 7 y2(0) 

In the model of fully ionized plasma the kinetic coefficients are also in a simple relation [5] 

(m) k 0.9 [m. 



n 



, ^«48.2>l; 
0.8 V m e 



(56) 



(57) 



the heat conductivity is determined by the fast thermal electrons while the viscosity by slowly moving protons. The 
ratio of kinetic coefficients is determined by the ratio of Maxwell velocities. 

By definition we have the boundary condition U(0) = 1 and for qualitative modeling we can suppose small 
derivatives of the wind velocity in the chromosphere region W(0) = 0. Analogously, we can suppose small gradient of 
the temperature F(0) = 0. Together with 0(0) = ©0 = s 2 /y~> 1 and T(0) = 1 we have the full set of boundary 
conditions at x = for the differential equations Eq. (l55t and Eq. (l44l for the wind variables U and ©. More 
complicated are the boundary conditions for the wave variables described in the next subsection. 



2.5. Boundary conditions for the waves 



With known background wind variables U (x) and T(x) we can solve the wave equation ( f20b for runaway AWs at 
x = I. As we will see later, the runaway boundary condition Eq. (l72l corresponds to right propagating AWs at the 
right boundary of the interval. The wave equation (l20l is extremely stiff at small viscosity, and a numerical solution is 
possible to be obtained only downstream from x = 0tox = l. We have to find the linear combination of left and right 
propagating waves at x = 0, which gives the runaway condition at x = I. 

The solution of wave equation according to Eq. (f23t determines the energy flux related to the propagation of AWs 



<7wave (¥(*)) 



I 

4 



-rjRe(u*w) 



where 



EL 

2no 



V -kn 



u 



Re 



u4 

2^0 





-In 
o 





(fa) 



(58) 



(59) 



Here the j'-term represent the kinetic energy of the wave, r\ -terms comes from the viscous part of the wave energy 
flux, and the Z?o-terms describe the Pointing vector of the wave. 

In order to take into account the boundary condition at x = I we calculate the eigenvectors of the matrix K, which 
according to Eq. (|20T > determine the wave propagation in a homogeneous fluid with amplitude °^exp(ikx). Then the 
eigenvalues of K give the complex wavevectors 



k = k' + ik" = eigenvalue(K), 



i.e. 



det(K-Jfcl) = 0. 

The three eigenvectors |L) , |D) and |R) are ordered by the spatial decrements of their eigenvalues 



k'l < < k'i < ftg, 



and are normalized by the conditions 



(Lk|L) 



(RlelR) = (Dip ID) 



1. 



(60) 
(61) 
(62) 

(63) 



where the sign corresponds to the direction of wave propagation. Notation L corresponds to left propagating wave, R 
to right propagating wave, and D for a overdamped at small viscosity mode. 
For technical purposes we introduce the matrix notations 



|L> 




|R> 




ID) 




(64) 



For low enough frequencies co — > and wind velocities the modes describes: 1) right-propagating AWs with kf R w 
0)/Va and small k R w Vk© 2 /2V^ ^ k' R , 2) left propagating wave = —k^, and a diffusion overdamped mode 
* V/t/VkU 3> k' D which describes the drag of a static perturbation by the slow wind U <C Va m a fluid with small 
viscosity. In this low-frequency and long-wavelength limit the stiffness ratio of the eigenvalues is very large 



N 
N 



vi 



k" 

k' ~ v k U co 



» 1. 



(65) 



The strong inequality is applicable to the chromosphere where the viscosity of the cold plasma is very low. As we 
emphasized, the wave equation (l20b is a very stiff system and indispensably has to be solved downstream from the 
chromosphere x = to the corona x = I using algorithms for stiff systems. Let 



!%»> = 




lv R M> = 




(66) 



are the solutions to the wave equation ( f20b with boundary conditions 

|VL(0)> = |L(0)), |V/ R (0)) = |R(0)). (67) 
We look for a solution as a linear combination 

V{x) = Vg.(x) + r VAx), (68) 

in other words we suppose that from the low-viscosity chromosphere plasma do not come overdamped diffusion 
modes. The strong decay rate make them negligible at x = 0. Physically this means that from the Sun are coming AWs 
(R-modes) and some of them are reflected from the transition layer (L-modes) 

Yr(0) = |R(0))+r|L(0)). (69) 

Analogously, for the configuration of open corona we have to take into account the runaway boundary condition for 
which we suppose zero amplitude for the wave coming from infinity 



Written by components 



iKZ)=f|R(Z))+c|D(/)>. (70) 



(71) 




this boundary condition gives a linear system of equation for the reflection coefficient r, transmission coefficient t and 
the mode-conversion coefficient c. 

For I — > °° when exp[— k^(l)l] <C 1 exp[—kf^(l)l] the amplitude of D-mode is negligible and the runaway boundary 
condition reads 

y(l) = y K (l) + ry L (l)*t\R(l)), (72) 

or by components 

*co W sxo W HIV (73) 



b R (i) J \b L (i) )~ v%(0 

These systems give the amplitudes of the reflected wave r and transmitted wave t in the solution Eq. d68l . For this 
solution we have the energy fluxes 

&=(Y(l)\g(l)\Y(l))= (74) 

= \t\ 2 + \c\ 2 + (tc*(D(l)\ g (l)\R(l))+c.c), 
l-^=( ¥ (0)\ g (0)\ ¥ (0))= (75) 

= l-M 2 + (r*(L(0)|g(0)|R(0))+c.c). 

Then we introduce the absorption coefficient 

(Y(x) \g(x)\ W (x)) \ l = \-M-3T. (76) 

The described solution is normalized by unite energy flux of the R-wave. If we wish to fix energy flux of the right 
propagating wave to be g W ave(0) we have to make the renormalization 

*F(x)=A wave lK*). (77) 
Now using we can calculate the wave part of the energy flux Eq. (l58l and the wave part of the momentum flux 

Yl v/ave (x) = -^-\b(x)\ 2 . (78) 
4-Hq 1 1 

This section is written in dimensional variables, but all equations can be easily converted in dimensionless variables 
as is done in the next sub-sub-section. 



2.5.1. Dimensionless wave variables 



Using mechanical units for length /, velocity f/o and density po we can convert all equations in dimensionless form. 
The formulas remain almost the same and we wish to mention only the differences. Introducing dimensionless density 



and wave energy flux 



we have dimensionless matrices 



p(x)=p(x)/p = l/I/(x) 

•^wave(O) = ~~m = e -^0 1^ wave | 5 



/ 

v{-iW+W) 

, mu -v\(-iw+W) (y 2 A -u 2 ^+u 2 djfi 



-vU 

— V 



and 



^4 



tt2 



1 

-Tj(x) 

,277^ 



-a 2 




2a 2 U(x) 











vi(o) = 


a 2 U(0) 



For the dimensionless energy flux we have 



1 i-,2 a 2 



U 



where 



U - 



Re(b g))--TfRe(§ #) = <^|f|^), 



^ Zd A M Zo 



Then for dimensional wave energy flux we have 

<7waveM — PqUq =2 waV e (x), 

and analogously for the momentum flux of the wave 

n wa ve(-^) = PqUq Sfi wave(-^); ^ wave 

in the next section we will consider all parts of the energy and momentum fluxes. 



(79) 
(80) 



(81) 
(82) 

(83) 



(84) 



(85) 



2.6. Mass, energy and momentum fluxes 



In the one-dimensional model which we now analyze, the conservation laws Eqs. (HJ-® are converted into three 
integrals of our dynamic system describing the mass j = poUoj, energy q = PqUq£?, and momentum II = PqUq^ 
fluxes 



j = Up = l, 



^U 2 +c p & Q T 



(86) 
(87) 



Jc&oF 
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- ( ^ 



-rf+£)UW 
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H 









b v -bv 



1 „ 

-r\ yv w + v 



w 



u 



const, 



_ 1 



-T7 + C )W + - 



const. 



(88) 



Here we can recognize the energy flux of an ideal inviscid gas 



q) ideal 
wind 



1 -U 2 + c p &qT, & = & Q T, 



(89) 



dissipative energy flux of the wind related to heat conductivity and viscosity 



^ s n s d = -x0 o F-(-77 + C)C/W, 



- /4 



(90) 



the non-absorptive part of wave energy flux 



6) ideal 
•^wave 
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2 -Re(b*%j\ 


4 











(91) 



and the absorptive part of the wave energy flux 



« e = --TfRe(iJW 



(92) 



proportional to viscosity. Analogously, for the momentum flux we have: 



£3 ideal _ 77 i ®® T 
^wind — ^ "T JJ ' 
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wind 



3 ideal 
wave 
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(93) 

(94) 

(95) 
(96) 



for the transversal AWs. As a rule, the dissipative fluxes are against the non-dissipative ones. One can introduce non- 
ideal wind energy flux 
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and non-ideal wind momentum flux 
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Then according to Eq. d28l l we have 



t(x) = 
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3. SHORTWAVELENGTH WKB APPROXIMATION 



It is instructive to apply short-wavelength approximation when we have weakly damped plane AWs. The secular 
equation (f2Tb 

V 2 A k 2 = C0d(c0d + iv k k 2 ) , (Oq = co-kU, (101) 

can be solved by successive iterations. In initial approximation Vk = for low frequencies we obtain for the right 
propagating wave (R-mode) k' m 0)/(Va + U). Then in linear with respect to Vk = approximation we derive 

2k w — — = ro<*- ( 1Q 2) 

y A V A (V A + U) 2 

Attenuation of the wave amplitude oce - *'* describes the damping of the wave energy density °ce~ 2k " x , or 

d.r£\vave = 2& £wave- (103) 

The wave energy flux obeys the same equation if additionally the dissipative wind energy flux is much smaller than 
the wave energy flux 

'<|^wave|, (104) 



^wind 



.e., if heat conductivity and wind viscosity are negligible. Formally this means that in Eq. d87l i F = = W, i.e. 



d x T 



d x U 



U 



which is the condition for the applicability of the WKB approximation — the medium should be approximately 
homogeneous. According to Eq. d98l l for U = const the momentum flux is the time averaged pressure of the AW 

n«V A = ^ = e fo , (105) 

which is simultaneously equal to the time averaged magnetic field energy density e b of the wave. For AW the wave 
part of the time averaged density of the kinetic energy is equal to the magnetic one 

e« « \p (« 2 ), = \p l«| 2 = £b- (106) 

For the averaged energy density of AW we have £ waV e = £» + £/>■ Then if Eq. d97l > we neglect all dissipative terms 
supposing homogeneous medium for a plane AW we obtain 

q{ X ) W (V A + C/)£wave + Ue u . (107) 

Let us summarize the physical conditions: we suppose runaway waves with small damping, there are no reflected 
waves by inhomogeneities. This result, notably Eq. (1107b . should be expected: the wave energy £ wa ve running with 
AW velocity V A is also drifting with the wind velocity U. This is an example of the aether drag theory, the energy 
density is moving with velocity V A + U. Additionally, the kinetic energy of the wave, e„, is blown away by the wind 
velocity U. 

Contrary to this, for the momentum flux there is no drag, since Maxwell tensor is not affected by the wind n « £/,. 
A simple integration of Eq. ( 1103b gives q(x) — qoD(x), where k"(x) is taken from Eq. (1102b 

D(x) ee exp (-2 rk"(x)dx) = expi - f <o 2 v k (x)dx 1 

V J F \ Jo [V A (x) + U(x)] 2 V A (x)j 

Rewriting Eq. ( 1107b as 

q{x) = (V A + C/)e W ave + Ue u = (2V A + 3U)s b , (109) 
and using Eq. ( 1105b we obtain the relation 

fls Sv rr^TT? v (110) 

2V A {x) + 3U(x) 



For the boundary conditions at x = we suppose that the wind velocity Uq is much smaller than Alfven velocity Va(0) 

J(0) = (2s + 3)#(0)> #(0), (111) 



i.e. in this dimensionless units the energy flux of AWs is much greater than momentum flux. Roughly speaking, in the 
beginning AWs only heat the corona rather than accelerate the wind. 

As main qualitative advantage of the MHD calculation of the solar wind problem is the natural explanation of the 
narrow width of the transition layer X . Formally we can define 



1 

I 



max 



d x T{x) 



T(x) 



(112) 



This width can be calculated using the WKB approximation for the energy and momentum fluxes. 



4. SELF-CONSISTENT PROCEDURE 

First we fix the boundary conditions, namely the temperature Tq and proton density n p (0) for x = 0. For these 
parameters we calculate the density po, Debye radius fb(0), Coulomb logarithm Jzfo, viscosity rjo, heat conductivity 
kq, Ohmic resistivity pn.o, and sound speed c s (0). The initial velocity of the wind is better to be parametrized by 
the dimensionless parameter s ^> 1, i.e., Uq = c s (0)/s. Analogously, plasma beta parameter j3o determines the Alfven 
speed at Va (0) = J y/2/3c s (0) . Let us also fix the maximal frequency CO for which we will consider plasma waves and 
calculate the absorption rate of the energy density of Alfven waves 2k" '(0). One can choose the interval of the solution 
to the MHD equations to be much larger than the AWs mean free path 1 /2k" (0), for example 



10 iov£(o) 



2k"(0) v k (0)co 2 



(113) 



Having units for length /, velocity Uq and density po, we can calculate dimensionless variables at x = : xq, tj , and 
0o- 

The input parameters of the program are Tq, « to t(0), Po, CO, and s which parametrizes j = (m) n tot (Q)Uo. We calculate 
I, a, Tj = Vo, Kq an d choose some A wave which finally determines increasing of the temperature T(l) = T(l)/T(0). 

The first step in the self-consistent procedure is to choose some initial approximations for the wind variables, the 
simplest possible functions are T(x) = U(x) = l and of course W(x) = F(x)= 0, i- e -> to extend the boundary conditions 
at x = for the whole interval x 6 (0, 1). Then one can start the successive approximations. 

1 . At fixed wind profiles T (x) and U (x) we calculate W , F , A = Ao + | In T + i In U , rf , Jc, d^ff, and dimensionless 
matrices M and g. Then we have to solve the wave equation, and to renormalize the solution with some fixed 
dimensionless energy flux for the R-mode J2 wave (0). 

2. Using so obtained wave variables *P we have to solve the ordinary differential equations for d^F and djW, with 
boundary conditions, for example, F(0) =0 and W(0) = 0. The variables 77 , Jc, djTj,, and djx, which participate 
in the coefficients have to be calculated simultaneously. 

3. Having solved the equations for the wind variables and formerly the equations for the wave variables we can 
calculate the total energy and momentum flux. If these fluxes are not constant with the required accuracy we go 
to item 1 and repeat the procedure. 

Now we can calculate the width of the transition layer by the maximal value of the logarithmic derivative of the 
temperature 

I , -, 

- =max d-lnT . (114) 

A x 

The width of the transition layer X, and the increasing of the temperature T(l) are functions of the wave energy flux 
coming from the chromosphere g W ave(0). We have to repeat the calculations for different A wave in order to derive the 
required temperature enhancement r(l;A wave ) = T{l)/T(0). At fixed enhancement and chosen spectral density of the 
AWs we have no more free parameters and have to compare the calculations with another model for heating the solar 
corona. 



5. CONCLUSIONS 

In spite that iron was suspicious from the very beginning, the problem of Coronium was a 70-year standing mystery 
until unambiguous identification as Fe 13+ by Grotrian and Edlen in 1939. The same 70-year time quantum was 
repeated. In 1947 Alfven J2[| advocated the idea that the absorption of AWs is the mechanism of heating the solar 
corona. Unfortunately the idea by Swedish iconoclast H was never realized in original form: what can be calculated, 
what is measured, what is explained and what is predicted. That is why there is a calamity of ideas still on the arena, 
for a contemporary review see the SOHO proceedings H. From a qualitative point of view the narrow width of the 
transition layer w = min |dx/dlnr(x)| is the main property which should be compared against the predictions of other 
scenarios. For example, in order the nanoflare hypothesis to be vindicated fH] such reconnections need to explain the 
narrow width of the transition layer at the same boundary conditions of wind velocity and temperature. Moreover, 
electric fields of the reconnections heat mainly the electron component of the plasma. How then proton temperature 
in the corona is higher? Launching the Hinode gave a lot of hints for the existence of AWs in the corona U, see also 
fiol [Til [l2ll . However most of the those research was in the UV region when high-frequency AWs which heated are 
already absorbed. All observations are for low-frequency (mHz range) AWs for which the hot corona is transparent. 
The best can be done is to extract the low-frequency behavior of the spectral density of AWs and to extrapolate to 
higher frequencies responsible for heating. So, observed AWs are irrelevant for heating. In order to identify AWs 
responsible for heating, it is necessary to investigate the high-frequency (1 Hz range) AWs in the cold chromosphere 
using optical, but not UV, spectral lines. The first hint in that direction is already found: 20 mHz AWs have been 
detected ifTHl . However, we are unaware whether new experiments of such a type are planned. One of the purposes 
of the present work is to focus the attention of experimentalists on the 1 Hz range AWs in the chromosphere, which 
we predict on the basis of our MHD analysis. For such purposes we suggest Doppler tomography lfl3ll of H a to be 
used; Ca lines are anot her p ossibility. Doppler tomography was successfully used for investigating rotating objects, 
such as accretion disks lll4ll and solar tornados lfl5l . Here we wish to mark also the Doppler tomography by Coronal 
Multichannel Polarimeter build by Tomczyk lfl6ll . The first observations of periodic variations of full-width-at-half- 
maximum of H a line profile also gave promising results lHoll . For investigating AWs by Doppler tomography we 
suggest the development of frequency dependent Doppler tomography operating as a lock-in voltmeter. The data from 
every space pixel should be multiplied by sincof and integrated for many wave periods. Finally, one can observe 
the time averaged distribution of the AW amplitude. Systematic investigation of such frequency dependent Doppler 
tomograms will reveal that Swedish iconoclast |0] is again right that the AWs heat the solar corona, after another 70 
years of dramatic cornucopia of ideas. 

ACKNOWLEDGMENTS 

One of the authors, T. M., is thankful to Yana Maneva for discussing many problems related to solar physics lfl7ll and 
Martin Stoev for the early stages of the present research lH8ll when the idea of self-induced opacity was advocated. 
He also appreciate fruitful discussions with Eckart Marsch on the problem of the high-frequency Alfven waves 
and hospitality of Sami Solanki at the MPI for Solar System Research. This work was partially supported by the 
St. Clement of Ohrid University at Sofia under Grant 203/2010. 

REFERENCES 

1. B. Edlen, Arkivf. Matem. Astr. och Fys. B 28, 1-4 (1941); B. Edlen, Zeits. f. Astrophysik 22, 30-64 (1942); P. Swings, 
Astmphys. J. 98, 116-128 (1943). 

2. H. Alfven, 0Mo/i. Not. R. Astron. Soc. 107, 211-219 (1947). 

3. H. Alfven, Nature 150, 405-406 (1942). 

4. R. P. Feynman, R. B. Leighton, and M. Sands, The Feynman Lectures on Physics: Exercises, Addison- Wesley, London, 1965, 
Problems 40.4-5. 



2 In this paper entitled "Granulation, magnetohydrodynamic waves and the heating of the solar corona" Alfven has suggested the heating mechanism 
of the solar corona. We note that this pionering work regretably is practically never cited in contemporary reviews on solar corona heating and 
lounching of the solar wind. 



5. L. D. Landau and E. M. Lifshitz, "Mean free path of plasma particles" in Physical Kinetics, 2nd edition, Butterworth- 
Heinemann, Oxford, 2002, sec. 43, Eq. (43.8-10). 

6. A. J. Dessler,B Science 170, 604-606 (1970); 

|http : //www . sciencemag ■ org/cgi/ content/ citation/ 17 0/3958/ 604] 

7. Proceedings of the SOHO 15 Workshop — Coronal Heating (6-9 September 2004, St. Andrews, Scotland), edited by 
R. W. Walsh, J. Ireland, D. Danesy, and B. Fleck, ESA SP-575, ESA Publications Division, Noordwijk, 2004. 

8. Ch. Day, Phys. Today 62, May 18-21 (2009). 

9. B. De Pontieu, S. W. Mcintosh, M. Carlsson, V. H. Hansteen, T. D. Tarbell, C. J. Schrijver, A. M. Title, R. A. Shine, 
S. Tsuneta, Y. Katsukawa, K. Ichimoto, Y. Suematsu, T. Shimizu, and S. Nagata, Science 318, 1574-1577 (2007); 
T. J. Okamoto, S. Tsuneta, T. E. Berger, K. Ichimoto, Y. Katsukawa, B. W. Lites, S. Nagata, K. Shibata, T. Shimizu, 
R. A. Shine, Y. Suematsu, T. D. Tarbell, and A. M. Title, Science 318, 1577-1580 (2007); R. Erdelyi and V. Fedun, 
Science 318, 1572-1574 (2007); Y. Katsukawa, T. E. Berger, K. Ichimoto, B. W. Lites, S. Nagata, T. Shimizu, R. A. Shine, 
Y. Suematsu, T. D. Tarbell, A. M. Title, and S. Tsuneta, Science 318, 1594-1597 (2007). 

10. D. B. Jess, M. Mathioudakis, R. Erdelyi, P. J. Crockett, F. P. Keenan, and D. J. Christian, Science 323, 1582-1585 (2009). 

11. J.-S. He, C.-Y. Tu, E. Marsch, L.-J. Guo, S. Yao, and H. Tian, Astron. Astrophys. 497, 525-535 (2009). 

12. B. N. Dwivedi and A. K. Srivastava, Current Sci. 98, 295-296 (2010). 

13. T. Marsh, Doppler Tomography, Europhysicsnews 36 No. 4, 133-138 (2005). 

14. T. R. Marsh and K. Horn, Mon. Not. R. Astron. Soc. 235, 269- 286 (1988). 

15. |http : / /www .mps .m pg . de/ en/ f orschung/ sonne/| 

16. S. Tomczyk and S. W. Mcintosh, Asprophys. J. 697, 1384-1391 (2009). 

17. Y. G. Maneva, "Generation and Dissipation of Alfven-Cyclotron Turbulence in the Solar Corona and Solar Wind and Related 
Ion Differential Heating and Acceleration," Ph.D. Thesis, Gottingen University, 2010. 

18. T. M. Mishonov, M. V. Stoev, and Y. G. Maneva, Eur. Phys. J. D 44, 533-536 (2007). 



3 Swedish Iconoclast Recognized after Many Years of Rejection and Obscurity 



